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We report evidence that empirical data show time varying multifractal properties. This is obtained 
by comparing empirical observations of the weighted generalised Hurst exponent (wGHE) with time 
series simulated via Multifractal Random Walk (MRW) by Bacry et al. [E.Bacry, J.Delour and 
J.Muzy, Phys.Rev.E 64 026103, 2001]. While dynamical wGHE computed on synthetic MRW series 
is consistent with a scenario where multifractality is constant over time, fluctuations in the dynamical 
wGHE observed in empirical data fail to be in agreement with a MRW with constant intermittency 
parameter. This is a strong argument to claim that observed variations of multifractality in financial 
time series are to be ascribed to a structural breakdown in the temporal covariance structure of stock 
returns series. As a consequence, multi fractal models with a constant intermittency parameter may 
not always be satisfactory in reproducing financial market behaviour. 
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I. INTRODUCTION 

The concept of multifractality in the context of fi- 
nance has received much attention in the econophysics 
literature particularly over the last two decades. Many 
empirical studies have investigated financial data scaling 
behaviour [IHZ] and several models have been proposed 
to account for the observed multifractal features [8Hl5|. 
Multifractal behaviour has consequently become a down- 
right stylised fact of financial market data [TH] , being ob- 
served across several classes of assets: from daily stock 
returns to foreign exchange rates and composite indices 
[2J [5] [TThT^] . If one denotes pt the asset price at time t 
and r tjT = pt+ T —pt the return at time t and scale r, then 
the process r tjT is multifractal if the following scaling is 
observed [5] 

E(|r tlT |«)=c s r«« ) (1) 

where c q is a constant and £ 9 is a non-linear function, 
called the scaling function. The departure from linear- 
ity of C, q is what distinguishes multifractal processes from 
uni-scaling ones. 

Multifractality is also particularly appealing for mod- 
elling financial markets as it offers a simple behavioural 
interpretation: looking at the volatility at different time 
scales is a very natural way to assess the impact of hetero- 
geneous agents in the market and therefore any measure 
of scaling behaviour can convey information about the 
efficiency of a given market [5D]. In this regard some 
authors have suggested and confirmed through extensive 
empirical studies, that scaling exponents can be repre- 
sentative of the stage of development of a market [2lll7j. 
The same works have convincingly shown that, through 
a hierarchy of scaling exponents, it is possible to classify 



markets according to their degree of development: emerg- 
ing (developed) markets exhibit scaling exponents signif- 
icantly larger (smaller) that those expected for purely 
uncorrelated noise. 

In the aforementioned spirit, the time evolution of the 
scaling function can be useful to track time varying prop- 
erties of the market. For this reason £ g has also been 
studied dynamically via the time dependent Hurst expo- 
nent [U [2H . In a recent publication [23] the authors 
have observed large fluctuations in multifractality mea- 
sured via the generalised Hurst exponent in empirical 
daily data [33] across different stock sectors. When one 
is faced with the task of ascertaining the nature of dy- 
namical fluctuations in these quantities, the subtle issue 
is being able to distinguish between bare statistical fluc- 
tuations, which are due to the finiteness of the sample, 
and true structural changes in the underlying multifractal 
process. In this paper we study the problem of validating 
dynamical fluctuations of the scaling functions, perform- 
ing an extensive analysis of of empirical stock returns and 
comparing their properties with those of synthetic multi- 
fractal series. Among all existing models we focus on the 
Multifractal Random Walk (MRW) introduced by Bacry 
et al. [8] because of its parsimonious formulation and its 
success in the econophysics literature. 
The paper is organised as follows: in Section [H] after 
giving a short summary of the main tools used in this 
context, namely the dynamical GHE method and an 
overview of the MRW model, we illustrate the relation- 
ship between the scaling exponents and the spectrum of 
the model. In Section |III| we report the main findings 
on the varying multifractality of empirical stock returns 
data, which are backed up by a simulation study on syn- 
thetic MRW series in Section [TV] a summary and conclu- 
sive remarks are drawn in Section fVl 
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II. GENERALIZED HURST EXPONENT AND 
MULTIFRACTAL RANDOM WALK 

The MRW model [5] is constructed, in its discrete ver- 
sion, via the sum of stationary increments 



{/At 

k=l 



(2) 



where each contribution is defined in the following way: 

x A t.k = e At , k e u, ^ k . (3) 

A continuous multifractal model is obtained in the limit 
At — > 0, but we will only consider the discrete case. In 
equation (|3j) eAt.k is a Gaussian white noise with vari- 
ance <7 • At and e" At fc represents the stochastic volatility, 
uncorrelated with e. By taking u)At,k as a stationary 
gaussian process, we have log-normal volatility compo- 
nents. What distinguishes the limit At — > of X\t{t) 
from a Brownian motion is the choice of the correlation 
structure of the process WAt,t, which is chosen, according 
to cascade-like processes [2"3] , as 



Cov(u;fc, Wfc+h) = A 2 log 



T 



{l + \h\)At 



(4) 



as long as it remains positive, and zero otherwise. The 
logarithmic decay with lag h of the auto covariance cre- 
ates long memory in the process. This specification im- 
plicitly defines the correlation length T and the inter- 
mittency parameter A. The appeal of this model for 
describing stock returns evolution lies in its ability to 
reproduce faithfully the most common stylised facts of 
financial markets. Its stationarity and its causal struc- 
ture make it also preferable to other multifractal models 
[25] . Among the three parameters of the model A, T and 
a, A is the one governing the non-linearity of the Cq spec- 
trum of MRW via the relation [5] : 



( q = (q-q(q-2)\ 2 )/2. 



(5) 



The following relation between the generalised Hurst ex- 
ponent (see Appendix [Alfor details) and the spectrum of 
£ g defined in equation (151) can be established 



C g =qH(q). 



(6) 



Since H(q) depends only on A in equation d6b s any time 
dependence for H (q) should be reflected in a time depen- 
dence for A. By considering q — 1,2,3,4, equations ([5| 
and Q the following expressions 



ff(l)=-(l + A 2 ), 
ff(3)=^(l-A 2 ), 



H{2) 
H{A) 



(1 + 2A 2 ). (7) 



Values of these expressions for different A's are reported 
in table [T] We perform the following Monte Carlo study 



TABLE I. Values of H(q) for q = {1,2,3,4} computed 
through the expressions in |7f . 





A = 0.2 


A = 0.3 


A = 0.35 


H(l) 


0.52 


0.5450 


0.5613 


H(2) 


0.50 


0.50 


0.50 


H(3) 


0.48 


0.4550 


0.4387 


H(4) 


0.46 


0.41 


0.3775 



in order to compare theoretical predictions with GHE's 
computed on synthetic time series: for a given value of A, 
we simulate 5000 MRW series and we compute the corre- 
sponding GHE's for every synthetic series; then, assum- 
ing these GHE's are i.i.d., which is quite reasonable since 
every MRW simulation is independent from the others, 
we obtain a sample distribution for the GHE's. From 
this distribution we compute the {2.5%, 50%, 97.5%}- 
quantiles which give the range of fluctuation of the GHE 
beyond which any observation can be deemed anomalous. 
We perform this study for q — {1,2,3,4}. Results of this 
analysis for three different values of A arc reported in 
table [Hj By comparing these values with those expected 
from the MRW model in table U we can conclude that 
all values computed via the equations given in ^ are 
compatible with the statistical confidence intervals com- 
puted on the synthetic time series. Nonetheless, one can 
see that confidence intervals for GHE are quite large and 
do not impose strict bounds to the values retrieved from 
the model (see a further discussion in Appendix [B]) . The 
quant iles for H{\) and H(2) reveal no relevant depen- 
dence on A, while a slight dependence is spotted on the 
2.5%-quantiles of H(3) and H(4), with the quantile de- 
creasing along with an increase of A. 



III. 



ANALYSIS OF EMPIRICAL DATA 



We now turn to show the main findings of the present 
paper, namely that stock returns exhibit time vary- 
ing multifractal properties. We analysed daily stock 
data quoted in the NYSE in the period ranging from 
01-01-1995 to 22-10-2012. The data were provided by 
Bloomberg. In [33] we argued that observed variations 
in GHE over time could be associated to a transition be- 
tween different regimes of stability for the firms studied. 
We have also considered as a measure of multifractal be- 
haviour the quantity 



AH w {q,q)=H w {q)-H w {q), q^q 



(8) 



where H w (q) is the weighted Generalized Hurst Expo- 
nent (wGHE) [33] . It was pointed out that the dynamical 
evolution of AH w (q,q ) (which we label AH™(q,q )) is 
useful to track possible changes in multifractality occur- 
ring over time. To test the significance of this statement 
with respect to a multifractal model for stock returns, we 
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TABLE II. We report numerical results of the quantiles obtained from the distributions of the GHE's computed on several 
synthetic MRW time series for H(q), q — 1,2, 3, 4. While -ff(l) and H(2) look quite stable across different values of A, one may 
note a dependence on A particularly in the 2.5%-quantile of H(3) and i?(4). 













{2.5%, 50%, 97.5%}quantiles 












£1(1 J 


H(2) 


H(3) 


£1(4 J 


A 


= 0.2 


T= 
T= 


=2500 
=5000 


{0.4256,0.5046,0.5870} 
{0.4293,0.5023,0.5749} 


{0.4260,0.4998,0.5794} 
{0.4197,0.4970,0.5736} 


{0.4127,0.4939,0.5799} 
{0.4107,0.4924,0.5726} 


{0.4050,0.4840,0.5852} 
{0.4039,0.4831,0.5762} 


A 


= 0.3 


T= 
T= 


=2500 
=5000 


{0.4281,0.5105,0.5911} 
{0.4266,0.5081,0.5839} 


{0.4152,0.4995,0.5867} 
{0.4241,0.4981,0.5892} 


{0.3936,0.4830,0.5803} 
{0.4017,0.4861,0.5806} 


{0.3831,0.4753,0.5891} 
{0.3778,0.4734,0.5722} 


A 


= 0.35 


T= 
T= 


=2500 
=5000 


{0.4271,0.5112,0.5975} 
{0.4212,0.4986,0.5830} 


{0.4167,0.4984,0.5918} 
{0.4280,0.5113,0.5856} 


{0.3891,0.4859,0.5935} 
{0.3899,0.4814,0.5801} 


{0.3653,0.4747,0.5832} 
{0.3594,0.4666,0.5684} 



defined in the introduction, and define 



- Empirical data 
-fBm with H=0.75 




y t = log|r t | = logcr t +log|e t | 



(9) 



FIG. 1. The scaling function for Microsoft Corp daily data is 
compared to that of a notorious unifractal model, i.e. the frac- 
tional Brownian motion. The choice of H = 0.75 means we 
are considering a process with long memory although unifrac- 
tal. 



proceed as follows: first of all we verify that the scaling 
functions of the data are indeed concave (see figure [l]for 
an example, where we report the scaling function of the 
daily stock returns of Microsoft Corp compared to that 
of a fractional Brownian motion). Then, after estimat- 
ing the parameters of the model from empirical data, we 
generate a large number of synthetic MRW time series 
and for each sample compute the GHE's for q = 1,2, 
from which we get AH W (1,2). This endows us with a 
sample distribution of independent AH W (1, 2)'s. From 
the cumulative distribution function of AH W (1, 2)'s one 
can extract confidence intervals. The latters provide the 
range of variation beyond which observed fluctuations 
can be considered to be significant with respect to the 
underlying model. 

The parameters of the model can be extracted from the 
empirical data through the behaviour of the volatility au- 
tocorrelation function. Consider r t .o — r t , where r t T was 



Then 



Vt - Vt+l = log 



( a t 



\ a t+i 

We can thus compute 



+ log|e t |-log|e t+i |. (10) 



((Vt - Vt+l) 2 ) = ( flog -^-) ) + ((log |e t | - log \e t+ i\) 2 ) 

(11) 

= C(l) + A(l-8 l>0 ). 

Therefore, ((yt — yt+i) 2 ) is proportional, up to a shift 
given by a constant A, to the log-volatility correlation 
function. In the MRW model we have for the log- 
volatility 



(12) 



and thus we can study ((yt — yt+i) 2 ) vs log/ to get an 
estimate of A 2 as well as of the correlation length T. A' is 
another constant which depends on the ratio T/At and 
therefore is used to estimate T. The other parameter of 
the MRW, i.e. a, can be computed from Var(r t ) = a. 
We thus estimate both T and A via linear fits of the 
((yt — yt+i) 2 ) as function of log/. Values of the param- 
eters obtained for different stocks are reported in table 
|III| together with the quantiles obtained from the em- 
pirical distributions of the GHE's. One can appreciate 
that the quantiles do not seem to vary much across dif- 
ferent stocks. In table III we also report the percentage 
of AH™ (1,2) exceeding the quantiles. The rather high 
number of points falling outside the confidence intervals 
confirm unassailably that empirical data systematically 
refuse the hypothesis of constant multifractality. In fig- 
ure [2] we show the results of this study for several daily 
stock prices, where AH™ (1,2) is shown to exceed the 
quantiles many times. Note that both IBM and Microsoft 
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TABLE III. The parameters of the MRW obtained for a number of stocks from our database, computed by means of the method 
described in the text. Note the very large values obtained for T, which are in principle difficult to justify given that the length 
of the time series is 3360 days. Neveretheless we have verified that shorter correlation lengths (up to the order ofW 2 ) do not 
make the extreme quantiles larger and hence do not jeopardise the validity of the results. 





A 


T (days) 


a 




Quantiles 


Exceedances percentage 


Boeing Corp 


0.1455 


608 


1.0651 


{- 


-0.0205,0.0039,0.0265} 


30 % 


Microsoft Corp 


0.2436 


739 


0.6064 


{- 


-0.0210,0.0072,0.0316} 


53 % 


PNC Financials 


0.2839 


538 


1.903 


{- 


-0.0206,0.0098,0.0349} 


15 % 


Sara Lee Corp 


0.222 


746 


0.2937 


{- 


-0.0204,0.0071,0.0318} 


47.7% 


Rowan Cos INC 


0.2046 


1080 


0.7992 


{- 


-0.0215,0.0066,0.0296} 


11.4% 


IBM Corp 


0.2363 


816 


1.7811 


{- 


-0.0214,0.0068,0.0310} 


39% 


Wells Fargo 


0.3490 


631 


0.5923 


{- 


-0.0237,0.0134,0.0433} 


9% 


American Express 


0.3278 


520 


0.8356 


{- 


-0.0212,0.0125,0.0429} 


14% 


General Motors Corp 


0.1727 


440 


1.0291 


{- 


-0.0211,0.0063,0.0278} 


41% 


Citigroup 


0.3798 


568 


6.4585 


{- 


-0.0202,0.0146,0.0458} 


13% 


JPMorgan Chase 


0.2950 


534 


0.9792 


{- 


-0.0219,0.0097,0.0356} 


27 % 


Bank of America 


0.3445 


670 


0.65 


{- 


-0.0256,0.0121,0.0420} 


30 % 


Morgan Stanley 


0.2528 


538 


1.1185 


{- 


-0.0221,0.0079,0.0333} 


14 % 



show a sudden decrease of multifractality corresponding 
to the 2007-2008 financial crisis, in apparent anticorre- 
altion with JP Morgan's A7J t u '(l,2). This finding is in 
agreement with what found in [23], where we observed 
that stocks from the financial sector exhibited increasing 
multifractality corresponding to the credit crisis. 



Empirical quantiles reported in table III and in figure 
and [2] have been computed from a Monte Carlo simula- 
tion with 1000 realisations of MRW's of 4000 time steps 
with parameters extracted from the empirical data. Thus 
our proxy of multifractality appears to defy the hypoth- 
esis of a constant multifractal behaviour, as in the setup 
of the MRW model. The analysis has been carried out on 
a large set of data from different sectors of the market. 
A large number of stocks show patterns similar to those 
reported in figure [2j 

We ascribe the inadequacy of the MRW to capture these 
observed empirical facts to the choice of a constant value 
of A. As we have already mentioned above, a time varying 
H(q) requires, through equation ([6]), a time varying A. It 
seems in fact reasonable enough to suppose the scaling 
properties of financial time series to vary over time be- 
cause of the reflection of complex and varying economic 
constraints. 

We must remark that the values obtained for the cor- 
relation lengths T are very large compared to the time 
series lengths, a feature which is nonetheless commonly 
observed in multi-timescale volatility models [15 . The 
observed long-range nature of correlations seems to call 
for the introduction of non-stationary models, which of 
course would present a number of further inconveniences. 
However, multifractality is expected to be proportional 
only to A 2 and thus the values of AHl" 2 we measured 



should not be affected by the true T of the series. In 
fact, by taking equation ([6| for two different values of q 
we have: 

C q Jqi-( q2 /q2 = H( qi )-H(q 2 ) (13) 
which gives after some rearrangement 



A 2 



(14) 



Therefore we expect our proxy of multifractality to in- 
crease along with an increase of A, regardless of the mea- 
sured T. 

Note also that the limit of weak covariance, i.e. for 
A — > 0, corresponds to a process whose behaviour is ap- 
proximately unifractal. As one can appreciate from the 



quantiles reported in table III the larger is A, the larger 
are the GHE quantiles either, which reflect the irregu- 
larities present in the data. Note though that there is 
no correlation between the number of quantile-exceeding 
GHE's and A. This fact tells us that a simple increase in A 
would not necessarily make empirical data fit the model. 
It seems more reasonable instead to suppose the covari- 
ance structure of the log- volatility not to be constant over 
time, reflecting some underlying market mechanism. 



A. The effect of t-Student residuals 

Since the hypothesis of the residuals in equation 
([3]) being Gaussian is at odds with the observed fat tails 
of the stock returns distributions, we have carried out 
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t 



FIG. 2. The black thick line represents the dynamical evolution of AH™ (1,2) (AH™ (1,2)) for four stocks from different 
sectors of the market: Microsoft Corporation, IBM CO, Boeing Corporation and JP Morgan. AH™ (1,2) is evaluated on 
moving overlapping time windows of 1250 days, with a shift of 100 days. The first time window includes data from 1995 to 
2000 and the corresponding AH™ (1,2) is plotted at the end of the period. The horizontal lines represent the 2.5% (red dashed 
line), 50% (green dot-dashed line) and 97.5% (magenta continuous line) quantiles extracted from the distribution of GHE 
estimated from 1000 realisations of MRW with parameters extracted from the data via the method described in the text. The 
quantile exceedances take place in all four cases with AH™ (1,2) always increasing on average. This suggests that in the period 
under study, the volatility may have exhibited a more complex behaviour than that predicted by constant multifractality models. 



the same analysis by taking the residuals to be t-Student 
distributed, that is 



P(e) 



1 



R r(|) ( e 2 + a 2y 



(15) 



where v is the number of degrees of freedom and the pa- 
rameter a is related to the variance a 2 via a 2 — a 2 / {u— 2). 
The student distribution is known to better fit stock re- 
turns tails, if one considers v = 4 or 5 [16]. The results 
obtained with this specification of the MRW model are 
reported in figure |4j We still retrieve a significant frac- 
tion of the observations overpassing the extreme quan- 
tiles, although this is reduced with respect to the case 
in which the residuals are Gaussian. We also observe 
a systematic shift upwards of the sample distributions 



of the AH W (1, 2)'s, with the median now around 0.03. 
This tells us that, for a fixed covariance structure, the 
thickness of the tails dramatically rebounds on the mul- 
tifractal properties of the simulated scries, which have 
now a higher degree of multifractality, a feature which 
had also been reported in [S]. 



IV. DYNAMICAL WGHE ON MRW TIME 
SERIES 



As a necessary test for the empirical analysis per- 
formed in the previous section, we now address the is- 
sue of studying the dynamical wGHE on synthetic time 
series generated via MRW. This model envisages a con- 
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FIG. 3. We show how the quantiles of the wGHE distribu- 
tion are affected by varying the correlation length T. We note 
no particular dependence, which in turn makes the test for 
varying multifractality depend only on A. 



stant value for the scaling exponents and multifractality 
for given intermittency parameter A and for each value of 
q, which means that possible variations of multifractality 
parameters for MRW times series are to be ascribed only 
to statistical fluctuations and are thus supposed not to 
cross the confidence intervals retrieved from their distri- 
butions. For this purpose, we test the wGHE method 
on a large set of MRW series to verify that the scal- 
ing exponents retrieved from these synthetic data are in 
agreement with the hypothesis of constant multifractal- 
ity, which corresponds to fixing the set of parameters of 
the model {A, T, a}. In figure [5] we give examples of dy- 
namical evolution H w (l) and AiJ™ 2 that one observes 
for MRW series. 

A complete list of the quantiles obtained for different 



values of A is given in Table IV together with the per- 
centage of Ai/™ 2 exceeding the confidence intervals. As 
also mentioned in the previous section, we have verified 
that the choice of T doesn't affect the quantile estima- 
tion: this is illustrated in figure |3j where we acknowledge 
that varying T, with A fixed, doesn't affect the distribu- 
tion of AHf 2 - For this reason we fix T — 500 and vary 
A. As one would expect, variations of these quantities 
occurring over time are well within the confidence inter- 
vals expected from the model. A comparison between 
the third column of table |IV| and the last right column 
of table III shows that, across different values of A, the 



number of A™ 2 exceeding the quantiles is much larger 
for empirical data than it is for synthetic MRW series. 
This analysis serves as a confirmation that the wGHE 
method returns values of multifractality and scaling ex- 
ponents which are consistent with the hypothesis of the 
MRW and thus confirms further the anomalous nature of 
the fluctuations observed on empirical data. 



V. CONCLUSIVE DISCUSSION 

The scope of this paper has been to measure and vali- 
date variations of multifractality on a set of stock returns. 
As a benchmark we have considered the MRW, a par- 
simonious multifractal model that is able to reproduce 
faithfully many features commonly observed in stock re- 
turns. Since the dynamical estimations of multifractal- 
ity are made on relatively small samples, to validate the 
observed fluctuations in the degree of multifractality as 
truly significant, we have first made an extensive study of 
the distribution of the generalised Hurst exponents H(q) 
for different values of q. We have shown that the fluctu- 
ations observed in empirical data are truly crossing the 
extreme quantiles and we have checked that the GHE 
of MRW series fluctuates instead well within the confi- 
dence intervals expected for a multifractal process with 
constant degree of multifractality. The confidence inter- 
vals have been computed on many realisations of MRWs 
simulated with a set of parameters obtained from differ- 
ent empirical data series. We ascribe the disagreement 
between model and data to the intermittency parame- 
ter being constant. We have shown that increasing the 
intermittency parameters makes the tails of the GHE dis- 
tribution computed on the synthetic series thicker, while 
we found no relevant dependence with T. Our analy- 
sis thus suggests that a varying intermittency parameter 
may be the correct guess towards the inclusion of the 
observed empirical facts into the MRW model. A possi- 
ble extension, which we are investigating at the present, 
is the possibility of identifying switching points between 
multifractality regimes, a topic which could be steered 
along the research stream of regime switching in dynam- 
ical correlation [26j [27] . The idea is to look for a covari- 
ance structure of the form 



Cov(uj k ,uj k+h ) =Q(k - k*)\llog 
+6(fc* - k)X 2 2 log 



Ti 



(l + |fe|)Ai 
{l + \h\)M 



, (16) 



where 0(x) denotes the Heaviside function. The sub- 
tlety of this approach is nevertheless the inclusion of a 
new parameter k* defining the switching point between 
different covariance regimes which, although stimulating 
for its possible economic meaning, seems to be tricky to 
be measured. 

The finding of very large correlation lengths is, in this 
framework, probably the best clue one has to detect non- 
stationarity in financial time series. The very large mem- 
ory present in the volatility is an intriguing concept which 
could also be related to the way markets evolves with 
history: the hypothesis of having a varying multi scal- 
ing behaviour might suggest that memory is itself a time 
varying process whose prolongation is to be ascribed to 
the growing complexity of modern societies which finan- 
cial markets reflect. In this regard, a more detailed study 
of the spectrum of the covariance structure of the volatil- 
ity on very large time series is likely to reveal some new 
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FIG. 4. The dynamical evolution of AH W (1,2) for the four companies IBM Corp, American Express, Boeing Corp and 
Bank of America. Now the quantiles have been obtained by considering a MRW with t-Student residuals. The distributions 
of the simulated returns have thicker tails than in the original Gaussian case, and this results in the quantiles of the sample 
distributions of AH W (1,2) being much larger. Nonetheless, the fluctuations observed on empirical data are still overpassing 
the quantiles significantly: 25% (IBM), 21% (Bank of America), 14% (American Express) and 18% (Boeing). We also note a 
significant shift upwards of the median of the sample distributions of AH™ (1,2) . This is a clue that multifractality is reflecting 
the tails of the synthetic returns distribution, although the covariance structure stays the same. 



insights (see [35J for a relevant discussion on this matter). 
The GHE method has been proved to be the most effi- 
cient among other scaling measures [4]; nonetheless, it 
would be interesting to investigate whether it is possi- 
ble to retrieve convincing evidence of variations in mul- 
tifractality by means of other methods. The study of the 
unconditional distribution of returns at all time scales 
could also be performed by means of parametric meth- 
ods; these have the drawback of requiring a specific model 
to be chosen though and, given the non-stationary be- 
haviour of financial time series, may be inappropriate, 
apart from very few practical situations. When one needs 
to carry out dynamic maximum likelihood estimation on 
the empirical data over the time windows at different 
time scales, the shortness of the time windows, especially 
at large time scales, results in the estimator not to con- 



verge and in the results to be hardly interpretable. 
Starting from the results of this work, our main future 
aim is the development of a realistic model for stock dy- 
namics which incorporate these empirical findings, pos- 
sibly backed by a plausible economic ground about what 
triggers multifractality to change over time. 



Appendix A: Generalized Hurst exponent method 

The generalized Hurst exponent H (q) is defined when 
the following scaling is observed: 

K q (r) oc T qH ^ , (Al) 
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FIG. 5. Left: Dynamical evolution (black thick line) of the wGHE H w (q = 1) computed on a realisation of a synthetic 
MRW time series. The horizontal lines represent the 2.5% (red dashed line), 50% (green dot-dashed line) and 97.5% (magenta 
continuous line) quantiles extracted from the distribution of GHE estimated from 1000 realisations of MRW with parameters 
A = 0.2, T=500, (7 = 1. Although statistical fluctuations are present they stay within the confidence intervals corresponding to 
the chosen set of parameters. Right: Dynamical evolution o/ Ai/™ 2 on another realisation of MRW with same parameters as in 
the left panel. As for the scaling exponent H w (q = 1), the fluctuations are compatible with a constant multifractality scenario. 
In both panels the shaded area represents the standard error respectively on the wGHE and on Aff™ 2 - The latter is computed 
as std(AHl 2 ) = y/var{H™{V)) + var{H™{2)), assuming that cov(H w {l), H w {2)) = 0. 



TABLE IV. Quantiles and number of AH™ 2 exceedances obtained from several simulations of MRW with different A. Note that 
the percentage of A_ff™ 2 exceeding the quantiles is negligible. This proves that the dynamic wGHE is faithfully reproducing the 
constant multifractality envisaged by the model. In all simulations we fixed T = 500 and a — 1, after verifying that this choice 
does not affect the results: in fact, all different values of T and a which we have tried return the same kind of percentages of 
exceedance. 



A 


Quantiles 


Exceedances percentage 


0.12 


{-0.0234,0.0033,0.0253} 


% 


0.14 


{-0.0227,0.0045,0.0258} 


1 % 


0.16 


{-0.0221,0.0048,0.0276} 


1 % 


0.18 


{-0.0223,0.0052,0.0272} 


1 % 


0.20 


{-0.0203,0.0056,0.0302} 


% 


0.22 


{-0.0244,0.0066,0.0310} 


% 


0.24 


{-0.0201,0.0076,0.0314} 


5 % 


0.26 


{-0.0211,0.0078,0.0330} 


1 % 


0.28 


{-0.0238,0.0081,0.0334} 


% 


0.30 


{-0.0219,0.0096,0.0357} 


1 % 



where K q (r) are the qth-order moment of the increment 
of a process X(t), defined as [29] 



K q (r) 



(\X(t + r)-X(t)n 
(\X(t)\ q ) 



(A2) 



r can vary between 1 and T max and (•) denotes the sample 
average over the time- window. Equation (Al) provides a 



way to compute directly the GHE for any value of q, as 
long as the scaling is observed. 



Financial returns time series exhibit often very large fluc- 
tuations, corresponding to periods of abnormal activity 
or large scale financial crises. These fluctuations may 
severely flaw any statistical estimator, when the latter is 
computed dynamically and the number of observations 
is therefore reduced; to track the dynamical evolution of 
the GHE, an exponential smoothing for the moments has 
been introduced in |23j , with corresponding weights cho- 
sen as the best trade-off between length of the time win- 
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dow and sensitivity to possible changes in scaling prop- 
erties of the time series. Details on the optimal weights 
can be found in [331 130] • Equation ( A2 1 is then replaced 
by 



find in fact 



(\S(t + T)-S(t)\*)„(t) 

(\S(t)\") w (t) 



(A3) 



where (-) w denotes weighted average. Analogously the 
weighted GHE (wGHE) H w (q) is computed by fitting 
the linear relation obtained inverting the scaling relation 
observed for the weighted moments: 



(A4) 



The exponent H w (q) embodies the scaling behaviour of 
the corresponding q-momcnt when the time series is ex- 
ponentially smoothed. 



Appendix B: Validity of the scaling 

A more detailed comparison between the relations Q 
and the the GHE's obtained from simulated series with 
different values of A has also been performed. We have 
simulate many synthetic MRW time series with A vary- 
ing in [0, ll and for each we have computed the GHE. 
In figure (pi) we plot the relations we obtain for H(q = 
1, 3) vs A. The black solid line in both figures is the best 
quadratic fit obtained from the measured points. We re- 
trieve, as expected from equations Q, quadratic trends 
for both H(l) and H(3) vs A, although the actual slope 
is not consistent with that predicted theoretically. We 



ff(l) = -(l + aA 2 ), 



ff(3) 



b\' 



(Bl) 



with a = 0.1348 and b — 0.2374 in the figure. Repeating 
the experiment many times we obtain a — 0.1475 ±0.011 
and b = 0.2502 ±0.017, still very far from the theoretical 
values a = b = 1. This apparent disagreement between 
the GHE method and the theoretical scaling function is 
probably due to the growing amount of noise present in 
the synthetic series as one increases A. Better fits are 
retrieved in fact if we shrink the upper bound of A: fits 
restricted to values of A < 0.3 work already much bet- 
ter giving a = 0.7162 ± 0.012 and b = 0.5512 ± 0.018. 
Note also that financial data exhibit, most of the time, 
intermittency parameters A < 0.3 (see also [31]). Some- 
where around A ~ 0.25 though, both fits start not to 
be distinguishable from a flat trend at H(l) = 1/2 and 
H(2>) = 1/2, as expected for H(2). We conclude that 
the noise present in the GHE doesn't allow to estimate 
the intermittency parameter reliably from the structure 
functions. 

Nonetheless this apparent flaw of the method doesn't af- 
fect the validity of the empirical results presented in sec- 
tion [TTTJ since we chose to estimate the intermittency pa- 
rameter from the volatility autocorrelation rather than 
from the scaling functions, although both methods being 
in principle valid. 
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FIG. 6. We plot the H(q) vs A for q — 1 (left) and q = 3 (right). The dots refer to the H(q)'s values obtained from the 
MRW- simulated time series with A ranging in [0,1]. The black solid lines are the best fits to eguations\Bl\ 
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